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Universal quantum computation can be achieved by simply performing single-qubit measurements 
on a highly entangled resource state. Resource states can arise from ground states of carefully de- 
signed two-body interacting Hamiltonians. This opens up an appealing possibility of creating them 
by cooling. The family of Affleck-Kennedy-Lieb-Tasaki (AKLT) states are the ground states of 
particularly simple Hamiltonians with high symmetry, and their potential use in quantum compu- 
tation gives rise to a new research direction. Expanding on our prior work [T.-C. Wei, I. Affleck, 
and R. Raussendorf, Phys. Rev. Lett. 106, 070501 (2011)], we give detailed analysis to explain 
why the spin-3/2 AKLT state on a two-dimensional honeycomb lattice is a universal resource for 
measurement-based quantum computation. Along the way, we also provide an alternative proof that 
the ID spin-1 AKLT state can be used to simulate arbitrary one-qubit unitary gates. Moreover, 
we connect the quantum computational universality of 2D random graph states to their percolation 
property and show that these states whose graphs are in the supercritical (i.e. percolated) phase 
are also universal resources for measurement-based quantum computation. 
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I. INTRODUCTION 

The rules of quantum mechanics appear to perform 
certain tasks much more efficiently than those of classical 
mechanics. The most celebrated example is the factor- 
ing of a large integer by Shor's quantum algorithm [T] 
that offers exponential speedup over existing classical al- 
gorithms. Quantum computers that implement generic 
quantum algorithms can take form in various computa- 
tional models, such as the standard circuit model [2J, adi- 
abatic quantum computer [3] 0] , and quantum walk [5] , 
all of which proceed via the important feature of quan- 
tum mechanics — the unitary evolution. 

A different but equally powerful framework is the 
measurement-based quantum computation [SHH]- A par- 
ticular computational model within this class is One- Way 
quantum computation [5] which we subsequently denote 
by MBQC. It proceeds by single-qubit measurements 
alone on a highly entangled initial resource state [8TllO|. 
For MBQC, resource states that allow universal quantum 
computation turns out to be very rare [T2J , but examples 
do exist [T3HTB] . The first identified universal resource 
state is the 2D cluster state on the square lattice [5J [T3J . 
It was also shown that 2D cluster states defined on reg- 
ular lattices, such as triangular, hexagonal and Kagome, 
are also universal resources [TT]. Cluster states and re- 
lated graph states can be created by the Ising interaction 
from unentangled states [13] and they have been created 
with cold atoms in optical lattices [T7]. However, they 
do not arise as unique ground states of two-body inter- 
acting Hamiltonians [18] . although they can be an ap- 
proximate unique ground state [19]. However, by going 
beyond qubit systems and by careful design of Hamilto- 



nians, a few quantum states have been found that are 
both unique ground states and universal for MBQC [2"0T - 
I23j . This opens up an alternative possibility of creating 
universal resource states by cooling the systems. 

Independently of the development on quantum compu- 
tation, Affleck, Kennedy, Lieb and Tasaki (AKLT) con- 
structed a family of states that were ground states of 
isotropic antiferromagnet-like Hamiltonians 24-215]. In 
any dimension, AKLT states are ground states of par- 
ticularly simple Hamiltonians which only have nearest- 
neighbor two-body interactions, are rotationally invari- 
ant in spin space and shares all spatial symmetries of 
the underlying lattice. In particular, AKLT provided an 
explicit example of a one-dimensional spin-1 chain that 
has a finite spectral gap above the ground state, sup- 
porting Haldane's conjecture on integer spin chains with 
spin rotation symmetry [27 . These valence-bond states 
turned out to be the first examples of matrix product 
states (MPS) [38J and projected entangled pairs states 
(PEPS) [13 H5]. The use of MPS and PEPS also gives 
rise to a new perspective on MQBC [2][TB]. In partic- 
ular, it was recently discovered that the one-dimensional 
spin-1 AKLT state [2UES] can serve as resources for re- 
stricted computations |141 130] . i.e., implementation of 
arbitrary one-qubit rotations. The discovery of the re- 
sourcefulness of AKLT states creates additional avenues 
for its experimental realization [31] ■ and has instilled 
novel concepts in MBQC, such as the renormalization 
group and the holographic principle [321 133], However, 
to achieve universal quantum computation within the 
measured-based architecture a two-dimensional structure 
is needed. 

In Ref. [3T], Cai et al. considered stacking up ID 



AKLT quasichains to form a 2D structure. Their con- 
struction showed that the resulting state, even though it 
is longer an AKLT state, can provide universal quantum 
computation. Later independently by us [34 and by Mi- 
ayke [35], it was shown that indeed the 2D AKLT state 
on the honeycomb lattice provides a universal resource 
for MQBC. Here, expanding on our prior work |34j . we 
provide an alternative proof that the ID spin-1 AKLT 
state can be used to simulate arbitrary one-qubit unitary 
gates, and generalize the method and give detailed anal- 
ysis to the proof that the spin-3/2 AKLT state on a two- 
dimensional honeycomb lattice is a universal resource for 
measurement-based quantum computation. We do this 
by showing that a 2D cluster state can be distilled by 
local operations. Along the way, we have connected the 
quantum computational universality of 2D random graph 
states to their percolation property. We note that ex- 
tension of our approach using Positive Operator Valued 
Measure (POVM) and percolation consideration to com- 
putational universality have been successfully applied to 
a deformed AKLT model in Ref. [36] . 

The structure of the present paper is as follows. In 
Sec. [II] we discuss how to locally convert a ID AKLT to 
a ID cluster state. In Sec. [TTT] we outline and illustrate 
the method of how to locally convert the 2D AKLT state 
to a random graph state. We then give the general proof 
in Sec |IV| In Sec. [V] we show the quantum computa- 
tional universality of these graph states is related to the 
percolation of the graph and show how to convert these 
graph states to a 2D cluster state on a square lattice. 
We support our assertion with Monte Carlo simulations 
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in Sec. [VI] and conclude in Sec. |VII| In the Appendices, 
we use a different approach to obtain the probability of 
getting any POVM outcomes. 
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FIG. 1. The ID AKLT state (a) and the ID cluster state (b). 



then made at every site that maps the state of the two 
virtual qubits to their symmetric subspace, which is then 
identified as the Hilbert space of spin-1 particle: 

P v = \S Z = 1)(00| + \S Z = -1)<11| + \S Z = 0)(i>+\, (1) 



|tf+> = — (|01) + |10», 



(2) 



where \S Z = ±1,0} are the three S — 1 angular momen- 
tum eigenstates: S Z \S Z ) = S Z \S Z ). For convenience we 
shall take the periodic boundary condition, so that the 
last site of the ID chain is actually connected to the first 
site of the chain. Open boundary condition can be dealt 
with by attaching qubits at the ends. The ID AKLT is 
therefore given by 
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which is the unique ground state of the following spin- 
isotropic Hamiltonian with a finite gap [531 US] : 
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II. ONE DIMENSION 

We begin by investigating the ID AKLT state and how 
it can be locally converted to a ID cluster state. By doing 
so, we have thus proved the equivalence of the capabil- 
ity to simulate one-qubit unitary gates for both types 
of states. Many of the methods developed in this section 
can be extended to the more interesting case of 2D AKLT 
state on the honeycomb lattice. 



where S v denotes the vector of the spin operator at site v. 
This AKLT model provided strong evidence in support 
of the Haldane's conjecture [27] , 

On the other hand, the ID cluster state \C\d) also 
can be understood similarly by projecting virtual entan- 
gled pairs to physical spins, known as projected entangled 
pairs states (PEPS) [29] , where the virtual entangled pair 
is replaced by \<p H )e = |00) + |01) + |10) - |11) and the 
local projection is given by P% = |0)(00| + |1)(11|, giving 
rise to 



A. ID spin-1 AKLT state and ID cluster state 



(5) 



The ID AKLT state [23] can be understood by us- 
ing the valence-bond-solid (VBS) picture, as illustrated 
in Fig. [T^,. (1) First one regards a spin-1 particle at 
each site as consisting of two virtual spin- 1/2 particles 
(qubits), each of which forms a singlet with the virtual 
qubit on the neighboring site: \<j>) e = 1 01) e — |10) e , where 
the normalization is omitted, |0) = | f) an( l |1) = I l) 
are eigenstates of Pauli a z , and e denotes the edge that 
links the two virtual qubits. (2) A local projection is 



However, for our purpose, it will be useful to define equiv- 
alently the cluster state as the common eigenstate of the 
following operators: 



Z v -\X v Z v+ i \C 



( 1D )) = |C( 115 ) 



(6) 



for all sites v, where v ± 1 are the two neighboring sites 
of v on the chain. Note that for convenience we denote 
the three Pauli matrices by X = a x , Y = a y and Z = a Zl 
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and use the two notations interchangeably. Moreover, the 
choice of or eigenvalue is arbitrary, as the re- 

sulting states are related by local unitary transformation. 
The ID cluster state can be used to simulate one-qubit 
unitary operation on one qubit and is the basic ingredient 
in MBQC 0. 

In fact, the ID AKLT state has been shown to be able 
to simulate one-qubit unitary operation as the ID clus- 
ter state [TH [501 E3 by explicitly constructing one-qubit 
universal gates. It has also been realized that the spin-1 
AKLT state can actually be converted, via local oper- 
ations, to the ID spin-1/2 cluster state with a random 
length [37 . In the following section, we provide an alter- 
native method for the reduction of the ID AKLT state to 
a ID cluster state. This method will then be generalized 
later for the reduction of the 2D AKLT state. 



B. Reducing ID AKLT state to a ID cluster state 

As spin-1 Hilbert space is of dimensionality three, in 
order to convert to dimensionality two of a qubit, a pro- 
jection or a generalized measurement is needed. In the 
mapping P v in Eq. there is a two-dimensional sub- 
space spanned by \S = 1,S Z = 1) and \S = 1,S Z = —1) 
or equivalently by the two virtual qubits 1 00) and |11). 
One can therefore consider 



F, 



\S, = l)(S z = l\ + \S t 



as a projection that preserves a two-dimensional sub- 
space, where we suppress the label 5 = 1. However, 
what happens if the projection is not successful and it 
ends up in the subspace orthogonal to that spanned by 
\S Z = 1) and \S Z = -1}? To solve this "leakage" prob- 
lem, one takes advantage of the rotation symmetry and 
adds two more projections: 

F x = {\S X = 1)(S X = 1| + \S X = -1)(S X = -l|)/\/2, (8) 

Fy = (\Sy = l)(Sy = l| + \Sy = -l){Sy = -1|)/A (9) 

and notice the completeness relation in the spin-1 Hilbert 
space: 
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S=l- 



(10) 



The above F's constitute the so-called generalized 
measurement or POVM, characterized by {F^F a }. Their 
physical meaning is to define a two-dimensional subspace 
and to specify a preferred quantization axis x, y or z. In 
principle, the POVM can be realized by a unitary trans- 
formation U jointly on a spin-1 state, denoted by |^>), 
and a meter state |0) m such that 



U\i>)\0) m = J2F a \i>)\a) 



(11) 



where for the meter states {a\a') — 5 a _ a >. A measure- 
ment on the meter state will result in a random outcome 
a, for which the spin state is projected to F a \ip) [2J. 



Claim. We shall show that after performing the gen- 
eralized measurement on all sites with {a v } denoting the 
measurement outcome the resulting state 
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is an "encoded" ID cluster state. 

In the following we shall make use of the equivalent 
representation of the AKLT state by the virtual qubits; 
see Eq. Q, e.g., \S Z = 1) = |00) and \S Z = -1) = |11), 
where the r.h.s are two-qubit states. In this regard, we 
can think of F operators in terms of two-qubit operators: 

F, = (|00}<00| + |ll)<ll|)/V2, (13a) 
F x = (\++)(+ + \ + \--)(--\)/V2, (13b) 
F y = + \-i,-i)(-i,-i\)/V2, (13c) 

where |±) satisfy cr x |±) = ±|±) and |±z) satisfy cr y \±i) — 
±|±i). Thus, in terms of these F's the post-measurement 



state ( 12 ) is simply given by 



(14) 



Naturally as with \S Z = ±1), there is also the corre- 
spondence between the other two 5 = 1 states and the 
two-qubit states in x and y bases: \S X = 1) = | + +), 



-l|)A/2 (7) \S X = -1) = | - ->, \S V = 1) = and \S y 



-1 



| — i, —i). The use of qubit enables us to take the advan- 
tage of the stabilizer formalism [35], even though its use 
is not essential. 

First, let us explain the meaning of "encoding" used 
in the claim. Suppose two neighboring sites u and v 
have the same outcome a = z; see Fig. [2^,. Then the 
two-dimensional subspaces at sites u and v are both 
spanned by 1 00) and |11). However, the singlet state 
between the two virtual qubits connecting u and v dic- 
tates that the appearance of the basis states are anti- 
corrclatcd. For example, |(00)„) at site u cannot coexist 
with |(00)„) at site v. There are only two possible basis 
states for the two sites u and v: |"0") = |(00) U (11)„) 
and |"1") = |(11) U (00)„). These two states "encode" a 
logical qubit {|"0"), |"1")}. In terms of spin-1 notation, 



they are \S Z = 1,5 2 



-l)u,v and \S Z 



T,5 2 
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showing the antiferromagnetic properties of the AKLT 
state. From these two states, logical Z and X opera- 
tors can be defined: Z = |"0")("0"| - |"1")("1"| and 
X = | "0"><"1"| + |<T')("0"| (and thus the Pauli opera- 
tor Y = —iZX can be determined). In the same manner, 
for k consecutive sites with the same outcome a, only one 
qubit is encoded by the k physical spins with quantiza- 
tion axis being in the a-direction. We shall refer to these 
sites collectively as a domain. On the other hand, for 
two neighboring sites having different outcome a u ^ a v , 
the four combination \S au = ±l,5 a „ = ±1) can appear 
and each site is effectively a qubit. 

The above analysis can be expressed in terms of the 
stabilizer formalism. In the example that neighboring u 
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POVM 
outcome 
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X 




stabilizer 
generator 
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AlA.i-\-\U y Uy 



a„ a v 



TABLE I. The dependence of stabilizers and encodings for 
the random graphs on the local POVM outcome in the case 
of ID AKLT state. \C\ denotes the total number of virtual 
qubits contained in a domain. In the first line, i — 1 .. 2\C\ — 1, 
and in the third line i — 1 .. 2\C\. One choice of the sign is 
Xi = 1 if the virtual qubit i is on the odd number of sites 
(relative to, e.g., the left end) in the domain and Xi = —1 
otherwise. 



and v share the same outcome a = z (i.e., the domain 
consists of two sites u and v) , for the two virtual qubits of 
site u (denoted by the 1 and 2) we have <8> cr^F UtZ — 

F UtZ . This means a[ 1] ® a l z ] |V(K») = WK})>- Sim- 
ilarly, for site v (with two virtual qubits labeled as 3 

[31 [31 

and 4) we have a z <g> a z \ip({a v })) = \i/j({a v })). How- 
ever, because of the singlet between 2 and 3, we have 
— er[ 2 ' (g) az \tp({a v })) = \ip({a v })}. The above three op- 
erators {erf 1 ' <X> a[ 2 \ a[ 3 ^ <g> a[ 4 \ — o^ 2 ' ® crl 3 '} are called the 
stabilizer generators and they define the logical qubit ba- 
sis states: |"0") = |(00)„(ll) tI ) and = |(ll) tt (00) B ), 
as one can verify that they are the common eigenstates of 
these operators with eigenvalue +1. The stabilizer gener- 
ators are effectively identity operators in the logical-qubit 
Hilbert space. To define the logical Z operator, there are 
many equivalent choices: e.g., Uz \ <Jz \ — c[ 3 ' and — cr[ 4 '. 
Any of them can be taken to one another by multipli- 
cation of some combination of the stabilizer generators. 
To complete the logical qubit operators, the X operator 
can be taken as X = <t x ^ , which flips |"0") to 

| "1"), and vice versa. Other outcomes can be dealt with 
in a similar way, and these are summarized in Table |T] 
Two important properties are that (i) each domain can 
contain more than one physical qubit and is only one 
logical qubit; (2) the qubit basis depends on the shared 
outcome of the POVM. 

We remark that even though a domain may contain 
two or more sites one can perform projective measure- 
ment on all but one site in the basis defined by {\S a = 
1) ± \S a = -1}}, where a is the label of the POVM out- 
come for the domain. The domain is then reduced to a 
single site but still preserves the same degree of entan- 
glement with its neighbors. 

To show that the post-POVM state is an (encoded) 
cluster state, let us illustrate with the example shown in 
Fig. [2]d. Let us label the three sites by u, v and w re- 
spectively. Suppose the POVM outcomes on these sites 
are a u = x, a v = z, and a w — x, respectively. First note 
that — a^a^ commutes with F u _ x and — a^a^ com- 

[21 [31 

mutes with F wx . Note also that —a x a x is a stabilizer 
operator of the singlet between 2 and 3, but it does not 
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FIG. 2. Illustration of (a) encoding and (b)& (c) stabilizer 
operator. 



commute with F v<z . Similarly, —(rif e!? is a stabilizer 
operator of the singlet between 4 and 5, but it does not 
commute with F vz , cither. However, if we multiply all 
the above operators, we obtain 



(15) 



Because commutes with F v>z , due to the identity 

<7 X ®(7 X (|00)<00| + |11)<11|) 

= (|11)(00| + |00)(11|) 

= (|00)(00| + |ll)<ll|)£T B ®(7 a!> (16) 

K v is thus a stabilizer operator for the post-POVM state. 

— [21 — 

In terms of logical Pauli operators Z u = a x , Z w = 
ct! 5 ', X v = crlfVI 4 ', we arrive at the stabilizer operator 
Ky = ZjjXyZw- This is the stabilizer operator defining 
a linear cluster state; see Eq. 

As a further illustration, let us consider the same three 
sites in Fig. |2jj but with a u = x, a v = z, and a w = y, 
i.e., the last site has a different outcome a w — y than 
the above example. Because of this, one now considers 
— Oj^CTjf instead of — ojjfioa^ and can show that the fol- 
lowing operator is a stabilizer generator: 



(17) 



Now, we use the logical operators Z u = er^ 2 ' , Z w = Oy^ , 
X v = a4 3 ' , and Z v = and we arrive at 



Ky — Z u (iX v Z v )Z w — ZyY v Z v 



(18) 



Although the stabilizer operator K v is not of the canoni- 
cal form of the cluster-state stabilizer Z u X v Z Wl they are 
related by local unitary transformation that leaves Z v 
invariant. 
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C. General proof of ID encoded cluster state 

The examples in the previous section prepare us for 
the general proof that the post-POVM state is an en- 
coded ID cluster state. Consider Fig. [2j;, in which there 
are three blocks labeled by U, V, and W, that may con- 
tain multiple sites having same POVM outcome, a u , a v , 
and a w , respectively. Let us label the last virtual qubit 
in block U by I, the first virtual qubit in block W by 
r, and the virtual qubits in block V by 1, 2, ... , 2k. Be- 
cause a v 7^ a u and a v ^ a w , we can separate the proof 
into two cases: (1) a u = a w , just as the first example 
(a u , a v , a w ) = (x, z, x) given in last section; (2) a u ^ a w , 
just as the second example (a u ,a v ,a w ) = (x,z,y) given 
in last section. The proof given below is a straightforward 
generalization of these examples. 

Case (1). Let us define a = a w — a u . For the edges 
connecting V to U and to W, consider the two operators: 
— CTa'cri 1 ' and — cri 2fc 'eri r '. Denote by a the label such that 
Xy = (^'j^aa^ is the logical X operator for the block V. 
For the edges connecting virtual qubits inside V, consider 

[21 [3] \2k — 11 

the operator: a a a a ..a a . The product of these three 
operators can be verified to be the stabilizer operator for 
the post-POVM state: 

K v = a^M ] Xv^ k] ^ ] - (19) 

As a / a„, either a a — cr a or a a — ±ia a a av and thus 
either a a a a = 1 or T^c^ , with ±cr a „ being a logical Z 

2 

for block V (and Z = 1 from contributions of virtual 
qubits 1 and 2r ). Using the encoded Z for block U and 

— r/i — r r i 

W, i.e., Z\j — ±a a and Zw = ±0a , we have 

K v = ±ZuX v Z w (20) 

is a stabilizer operator. (The choice of ± depends on the 
convention; see Table [i]) 

Case (2). For the edges connecting V to U and to 
W, consider the two operators: —(Ta u Oai and — <7a„ cr^j,. 
Denote by a the label such that X v = ®fL x a]a is the 
logical X operator for the block V . For the edges con- 
necting virtual qubits inside V, consider the operator: 
Ca'ca' ..fTQ fc_1 ' . The product of these three operators is 
the stabilizer operator for the post-POVM state: 

Xvs^M^P^P^W. (21) 

As a u ytz a w , a is either equal to a u or a w and hence 
the product of cr^J a$ aa k ^ Ocffl becomes either trl^crl 1 ' or 

\2k] \2k] — 

<Ta (Ta w . Either of them is a logical ±iZy- Thus, we 
have 

K v = ±iZ u(Z v X v )Z w — ±ZuY v Z w (22) 

is a stabilizer operator. This concludes the proof that 
the post-POVM state is an encoded ID cluster state. 



D. Probability of a POVM outcome 

Given a set of POVM outcome {a v }, what is the prob- 
ability p({a v }) that this occurs? This is can be obtained 
from the norm square of the resulting un-normalized 
post-POVM state \tp({a v })), namely, 

P(W) = (^({«,})^({^})>/(4 1 StI4 1 St)- (23) 

Let us denote by |V| the total number of domains, which 
is the number of logical qubits and \£\ the total num- 
ber of edges connecting domains. As we consider the 
periodic boundary condition, trivially \£\ = \V\, ex- 
cept when all sites have the same POVM outcome, i.e., 
a v = a for all v. Note that this latter case can never 
occur if the total number of the original spins is odd, as 
the frustrated configurations | + 1, —1, +1, —1, +1) and 
| — 1, +1, —1, +1, —1) (with the first and last sites be- 
ing connected next to each other) cannot appear [5S] . 
It turns out that, barring the exception of zero proba- 
bility, p({a v }) ~ 2' v ''~' £ '. This is because for contract- 
ing (tp({a v })\ip({a v })) to compute the norm we need to 
evaluate \{a, /3|(|01) — |10))| 2 , where a, (3 can be any of 
the six possibilities: {0, 1, +, — , +i, —i}. The ratio of the 
above expressions in the case were a and /? belong to dif- 
ferent bases to the case where they belong to the same 
basis (thus a — is 1/2. In total, there are 2' y l 

terms of equal contribution to the norm square, each re- 
duced by a factor 2~I £ L This results in the probability 
p({a4)~2l^H £ l. 

For the total number of sites n being even, all the 3™ 
possible POVM outcomes can occur, each with proba- 
bility po, except for the three configurations (with all 
a v being the same) having probability 2po. Solving 
(3" - 3)p +3- 2p = 1, we obtain p = 1/(3" + 3). For n 
being odd, the three configurations with all a v being the 
same cannot occur. All other configurations occur with 
a probability 1/(3" — 3) each. For large n, it is a very 
good approximation to regard all configurations {a v } as 
occurring with equal probability and hence the resulting 
ID cluster state contains on average 2n/3 qubits, which 
agrees with the result in Ref. [37]. 



III. REDUCTION OF THE 2D AKLT STATE 

Now that we have understood the ID case, to show 
that the 2D AKLT state is a universal resource for quan- 
tum computation, we proceed in three steps. First, we 
show that it can be mapped to a random planar graph 
state \G(A)) by local generalized measurement, with the 
graph G(A) depending on the set A of measurement out- 
comes on all sites. Second, we show that the compu- 
tational universality of a typical resulting graph state 
\G(A)) hinges solely on the connectivity of G(A), and 
is thus a percolation problem. Third, we demonstrate 
through Monte Carlo simulation that the typical graphs 
G{A) are indeed deep in the connected phase. We remark 
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that extension of our approach using POVM and per- 
colation consideration have been applied to a deformed 
AKLT model in Ref. [15]. 

The AKLT state [Ml US] on the honeycomb lattice L 
has one spin-3/2 per site of C. The state space of each 
spin 3/2 can be viewed as the symmetric subspace of 
three virtual spin-l/2's, i.e., qubits. In terms of these 
virtual qubits, the AKLT state on C is 

|$aklt)= (g) Ps,v ® ^> e > ( 24 ) 

veV(C) eeE(C) 

where V(C) and E(C) denote the set of vertices and edges 
of C, respectively. Ps,v is the projection onto the sym- 
metric (equivalently, spin 3/2) subspace at site v of C: 

P s = |000)(000| + | 111) <111| + \W)(W\ + \W)(W\, (25) 

where 

|W) = ~(|001> + |010> + |100)), (26) 

|W) = ~(|110> + |101> + |011)). (27) 

The mapping between three virtual qubits and spin-3/2 
is given by: |000) O {3/2,3/2), |111) o |3/2, —3/2), 
\W) O |3/2,l/2) and \W) O |3/2, — 1/2). For an edge 
e = (v,w), \4>) e denotes a singlet state, with one spin 
1/2 at vertex v and the other at w. For illustration, 
see Fig. [3^,. The AKLT state is the ground state of the 
following Hamiltonian: 

Haklt = X/ ' + 243 ^ ' ^ 3 ^ + 243 ^ i ' ^ ^ ' 

edge (ij) 

where an irrelevant constant term has been dropped. 

Next, we give the definition of a graph state [39], to 
which we shall prove that the AKLT state can be locally 
converted. A graph state \G) is a stabilizer state [35] 
with one qubit per vertex of the graph G. It is the unique 
eigenstate of a set of commuting operators [15] . usually 
called the stabilizer generators, 

X v (g) Z u \G) = \G), Vv e V(G), (29) 

u£nb(v) 

where nb(u) denotes the neighbors of vertex v, and 
X = <j x , Y = <jy and Z = a z are the three Pauli ma- 
trices. A cluster state is a special case of graph states, 
with the underlying graph being a regular lattice; see 
e.g. Fig. [3Jd for the illustration. Any 2D cluster state 
is a universal resource for measurement-based quantum 
computation [SI ITT]. 

To show that the 2D AKLT state of four-level spin-3/2 
particles can be converted to a graph state of two-level 
qubits, we need to preserve a local two-dimensional struc- 
ture at each site. This is achieved by a local generalized 




FIG. 3. Illustrations of the AKLT state on the honeycomb 
lattice, a graph state and the 2D cluster state on a square lat- 
tice, (a) AKLT state. Spin singlets of two virtual spins 1/2 
are located on the edges of the honeycomb lattice. A projec- 
tion Ps, v at each lattice site v onto the symmetric subspace of 
three virtual spins creates the AKLT state, (b) A graph state. 
One qubit (i.e. spin 1/2) resides at vertex of the graph. One 
stabilizer generator of the form X v ®ue»b(i;) ^ u * s sri ° wn - ( c ) 
2D Cluster state is a special case of graph states, where the 
graph is a two-dimensional square lattice. 



measurement [2] , also called positive-operator- value mea- 
sure (POVM), on every site v on the honeycomb lattice 
C. The POVM consists of three rank- two elements 

F v , z = j|(|000)<000| + |111>(1H|) (30a) 

Fv, x = yf (| + ++><+ + +| + |- -){ 1) (30b) 

F v , y = + I — — — < — — «, — i I ) , (30c) 

which extend those in Eq. (JT3J) to three virtual qubits. 
Note that |0/1), |±) = (|0> ± \l))/y/2 and | ± i) = 
(|0) ± i\l))/V2 are eigenstates of Pauli operators Z, X 
and Y, respectively. Physically, F v>a is proportional to a 
projector onto the two-dimensional subspace spanned by 
the S a = ±3/2 states, i.e., \S a = ±3/2) {S a = ±3/2|. We 
have simply used the three-virtual-qubit representation, 
and it will be useful for our proof. The above POVM 
elements obey the relation Y,ve{x,y,z} H,uFv,u = Ps,v, 
i.e., project onto the symmetric subspace of three qubits, 
equivalently, the identity in S = 3/2 Hilbcrt space, as 
required. The outcome of the POVM at any site v is 
random, x, y or z, and it can be correlated with the out- 
comes at other sites due to correlations in the AKLT 
state. As we demonstrate below, the resulting quan- 
tum state, dependent on the random POVM outcomes 



A={a v ,veV(£)}, 

l*M)> = (g) I^aklt) = (g) F„,«. (g) |0) e 

tiEV(£) ueV(-C) e£E(C) 

(31) 

is equivalent under local unitary transformations to an 
encoded graph state \G(A)). The graph G(A) determines 
the corresponding graph state, and we show that it is con- 
structed from the honeycomb lattice graph by applying 
the following two rules, given A: 

Rl (Edge contraction): Contract all edges e G E(C) 
that connect sites with the same POVM outcome. 

R2 (Mod 2 edge deletion): In the resultant multi- 
graph, delete all edges of even multiplicity and con- 
vert all edges of odd multiplicity into conventional 
edges of multiplicity 1. 

These two rules are illustrated in Fig. |4j A set of sites in 
£ that is contracted into a single vertex of G(A) by the 
above rule Rl is called a domain, which we have already 
encountered in the reduction of ID AKLT state. Each 
domain supports a single encoded qubit. The stabilizer 
generators and the encoded operators for the resulting 
codes are summarized in Table [TT1 Below we demonstrate 
the post-POVM state |W(^4)) is a graph state and justify 
rules Rl and R2 with simple examples. 

Rule 1: Merging of sites. Physically, this rule derives 
from the antiferromagnetic property of the AKLT state: 
neighboring spin-3/2 particles must not have the same 
S a = 3/2 (or -3/2) configuration pMH2"rJ] . Hence, after 
the projection onto S a = ±3/2 subspace by the POVM, 
the configurations for all sites inside a domain can only 
be |3/2, -3/2, . . . ) or | - 3/2, 3/2, . . . ), and intuitively, 
these form the basis of a single qubit. This encoding 
of a qubit can also be understood in terms of the stabi- 
lizer. Consider the case where two neighboring POVMs 
yield the same outcome, say z; see Fig. [5^. As a result 
of the projections F UjZ and F VtZ (with u = {1,2,3} and 
v = {4, 5, 6} each containing three virtual qubits), the op- 
erators ZiZ 2 , Z 2 Zs, and Z^Z 5 , Z 5 Z 6 become stabilizer 
generators of the post-POVM state ^(.A)). In addition, 
the stabilizer — Z3Z4 of the singlet state | </>) 34 commutes 
with the projection F UtZ (&F VtZ , and thus remains a stabi- 
lizer element for |W(^4)). In summary, the stabilizer gen- 
erators are Z1Z2, Z2Z^, —Z3Z4, Z^Z^, ZqZq, giving rise 
to a single encoded qubit 

aKoooum^ + ^Cmuooo),,), 

which is supported by the two sites u and v jointly. 
We observe here the antiferromagnetic ordering [24H26] 
among groups of three virtual qubits. To reduce the sup- 
port of this logical qubit to an individual site of £, a mea- 
surement in the basis {|(000)„) ± |(111)„)} is performed. 
The resulting state is a|(000)„)±/3|(lll) u ), with the sign 
"±" known from the measurement outcome. This is the 
proper encoding for a domain consisting of a single site. 
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(a) Rule R1: 




FIG. 4. Graphical rules for transformation of the lattice C 
into the graph G(A), depending on the POVM outcomes A. 

(a) Graph Rule 1: A single edge between neighboring sites 
with the same POVM outcome i £ {x, y, z} is contracted, (b) 
Rule 2: Pairs of edges between the same vertices are deleted. 

(b) An example for the application of the rules 1 and 2 to 
a small honeycomb lattice. The alphabets (j 7^ k) inside 
the circle indicate the POVM outcomes, (c) An example to 
illustrate the applications of the graph rules. 

(a) (b) 



u 




V 



FIG. 5. Illustrations of the encoding and the cluster graph 
stabilizer, (a) If two neighboring sites u and v share the same 
POVM outcome, e.g., a u — a v — z, then collectively these two 
sites form a logical qubit. More generally, if a set of connected 
sites share the same POVM outcome, then these sites effective 
encode a logical qubit. (b) An example of four sites u, v, w, 
and c with POVM outcomes a c = z, and 

used to illustrate the stabilizer generator K c — ±X C Z U Z V Z W . 

Domains of more than two sites are thereby reduced to 
a single site in the same manner. 

To see that the state ^(.A)) is indeed equivalent under 
local unitary transformations to the encoded graph state 
\G{A)), we consider the example of four domains c, u, v, 
w, each consisting of a single site of C, where the POVM 
outcome is z on the central domain c and x on all lateral 
domains u, v and w; see Fig. [5}d. By similar arguments 
as above, the operator O = —XiXyX2X2'X^Xy is in 



the stabilizer of |\&(./4)). Using the encoding in Table [h] 
i.e., with the encoded Pauli operators X c = Z\Z%Z^, 
Z u = ±Xy , Z v = ±X 2 > , and Z w — ±Xy, we find that 
O = ±X C Z U Z V Z W which is (up to a possible sign of 
convention) one of the stabilizer generators defining the 
graph state. Intuitively, we see that each edge from an 
outer domain contributes to an encoded Z from that do- 
main and the operator restricted at the center domain is 
clearly not an identity nor an encoded Z as its POVM 
outcome differs from the outer ones. This gives rise to a 
stabilizer generator local-unitarily equivalent to the one 
given in Eq. ( 29 ) . 

Rule 2: Mod 2 edge deletion. By the above construc- 
tion, if two domains u, v are connected by an edge of 
multiplicity m, the inferred graph state stabilizer gen- 
erators will contain factors of X U Z V or X V Z U . We 
observe that Z 2 — I, from which Rule 2 follows. 

Generalizing the above ideas, it is straightforward to 
rigorously prove that for any POVM outcomes for any A, 
the state |^(^4)) is local equivalent to an encoded graph 
state \G(A)); see below. We shall denote by \G(A)) the 
graph state after reducing multiple sites in every domain 
to a single site, i.e., to the proper qubit encoding by 
domains, as the graph remains the same. 



IV. FROM THE AKLT STATE TO GRAPH 
STATES: GENERAL PROOF 

Let us recall the POVM to be performed on all sites 
can be rewritten as: 



p 1 — 

± v,z — 

^ v,x — 

TP — 

r v ,y — 



2 l\2 + Z\Z 2 I23 + Z2Z3 

3 2 2 ' 

2 J12 + X1X2 I23 + X2X3 

3 2~ ~~2 ' 

2 I12 + YiY 2 I23 + V2Y3 

3 2 2 ' 



(32) 



It turns out that, for any A, the state ^(.A)) is local 
equivalent to an encoded graph state \G(A)), with the 
graph G(A) constructed as follows. An edge (v,w) £ 
E{C) is called internal iff at the sites v and w the local 
POVM has resulted in the same outcome. The graph 
G{A) is obtained from the lattice graph C by (1) con- 
tracting all internal edges, and, in the resultant multi- 
graph, (2a) deleting all edges of even multiplicity and 
(2b) converting all edges of odd multiplicity into conven- 
tional edges of multiplicity 1. See Fig. 2 for illustration. 

In step (1) of the above procedure, several sites of C 
are merged into a single composite object C £ V(G(A)). 
Each such C is both a vertex in the graph G(A) and 
a connected set of same- type sites of C, i.e., a domain. 
Physically, in a domain of type a, we have antiferromag- 
netic order along the ±a-direction, because two neigh- 
boring spins never have the same S a = 3/2 (or -3/2) 
in the AKLT state [23]. The state of the domain con- 
tains only two configurations w.r.t. the quantization axis 



POVM 

outcome 


z 


X 


V 


stabilizer 
generator 
X 

Z 


Xi\i+iZiZi+i, 

XiZi 


XiXi+iXiXi+i 

<g> -H 

XiXi 


Xi Xi+iYi Yi+i 

<8> -S z t 



TABLE II. The dependence of stabilizers and encodings for 
the random graphs on the local POVM outcome. \C\ denotes 
the total number of virtual qubits contained in a domain. In 
the first line, i = 1 .. 3\C\ — 1, and in the third line i — 1 .. 3|C|. 
The honeycomb lattice C is bicolorable or bipartite and all 
sites can be divided into either A or B sublattice, V(C) = 
A U B. Then, one choice of the sign is Xi — 1 if the virtual 
qubit i G v G A and A; = —1 if i £ v' £ B. 



a: I + 3/2, -3/2, +3/2, . . . ) and | - 3/2, 3/2, -3/2, . . . ). 
Thus, it is effectively one qubit. 

The outlined construction leads to one of our main 
results: 

Theorem 1 For any A that specifies all outcomes 0/ 
POVMs on C, quantum computation by local spin-3/2 
measurements on the state \^{A)) can efficiently simu- 
late quantum computation by local spin- 1/2 measurement 
on the graph state \G(A)). 

Thus, the computational power of the AKLT state, as 
harnessed by the POVMs Eq. ( |32| , hinges on the connec- 
tivity properties of G(A). If, for typical sets A of POVM 
outcomes, the graph state \G(A)) is computationally uni- 
versal then so is the AKLT state. 

The proof proceeds in three steps. First we show that 
every domain C £ V(G(A)) gives rise to one encoded 
qubit. Second, we show that |^(A)) is, up to local en- 
coded u nitaries, equivalent to the encoded graph state 
\G(A)). Third, we show that the encoding can be unrav- 
eled by local spin-3/2 measurements. 

Step 1: Encoding. Consider a domain C C V(C). 
That is, on all sites v £ C the same POVM outcome 
a £ {x,y, z} was obtained. C contains 3|C| qubits. The 
projections F v ^ a on all v £ C enforce 2\C\ stabilizer gen- 
erators, c.f. Eq. ( [30| . Furthermore, choose a tree T 
among the set of edges with both endpoints in the domain 
C. Each edge (u,v) £ T contributes a stabilizer genera- 
tor — (T^o-(^ to the product of Bell states <8>e6_e(£) \4>)e- 



These stabilizers commute with the local POVMs (32) 
and therefore are also stabilizer generators for ^(Aj). 
Since |7~| = |C| — 1, in total there are 3|C| — I stabilizer 
generators with support only in C, acting on 3|C| qubits. 
They give rise to one encoded qubit. 

While the stabilizer generators for our code follow from 
the construction, there is freedom in choosing the en- 
coded Pauli operators. Table |TT] shows one such choice of 
encoding. 

Step 2: We show that |^(A)) is an encoded graph 
state. Consider a central vertex C c £ V(G(A)) and all its 
neighboring vertices C p £ V{G(A)). Denote the POVM 
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outcome for all £-sites v £ C C ,C M by a c and an, respec- 
tively. Denote by the set of £-edges that run be- 
tween C c and C M . Denote by E c the set of £-edges in- 
ternal to C c . Denote by C c the set of all qubits in C c , 
and by C M the set of all qubits in C M . (Recall that there 
are 3 qubit locations per £-vertex v £ C C ,C M .) We first 
consider the stabilizer of the state ® e eB(£) \4>)e- F° r 
any /j, and any edge e £ En, let u{e) [v{e)\ be the end- 
point of e in C M [C c ). Then, for all fj, and all e £ E^ 

the Pauli operators — ai u J" e ^ are in the stabilizer of 

®e£E(c) l^e- Choose b £ {x,y, z} such that b ^ a c , 
and let, for any edge e' £ E c , Vi(e'),v 2 (e') £ C c be 
qubit locations such that e' = (vi(e'),V2(e')). Then, 



for all e' £ E c , —a 



(•ui(e')) (v 2 (e)) 



is in the stabilizer of 



®eeE{C) 

tors, 



Therefore, the product of all these opera- 



Thus, finally, for all C c £ V(G(A)), 



Ca £ { ±R Cc 



with R = X, Y 



C„£V(G(A)) 



This is, up to conjugation by one of the local encoded 
gates Ic e , Zc c , exp ( ± in /4 Zc c ) , a stabilizer generator for 

the encoded graph state \G(A)). The code stabilizers in 
Table [Tl] and the stabilizer operators in Eq. (34) together 
define the state |^(^4)) uniquely. |^(^4)) is, up to the 
action of local encoded phase gates, an encoded graph 
state \G(A)). 

Step 3: Decoding of the code. We show that any do- 
main C £ V(G(A)) can be reduced to a single elementary 
site w £ V(£) by local measurement on all other sites 
» £ C, » / in. For any such v, choose the measurement 
basis B a , a £ {x, y, z}, as follows 



Or =± 



(u(e)) (v(e)) 
w a u u a u 



Me'))(v 2 (e'))\ 
°b °b 



is also in the stabilizer of 



(33) 



e£E(C) 



We now show that Oc c commutes with the local 
POVMs and is therefore also in the stabilizer of |^(^4)). 
First, consider the central domain C c . The operator Oc a 
acts non-trivially on every qubit in C c , Oc a \i =/= h for 
all qubits I £ C c . Furthermore, for all qubits I £ C c , 
Oc a \l 7^ Oal- Namely, if I £ C c is connected by an edge 
e £ E^ to for some fi, then Cc \i = Oal ^ <r2} 
(for all [i, a M ^ a c by construction of G(A)). Or, if 
/ £ C c is the endpoint of an internal edge e' £ E c then 
Oc a \i — ctj 7^ (Ta} (a c ^ b by above choice). There- 
fore, for any i,j £ C c , Oq c anticommutes with Oal and 
o~a}, and thus commutes with all Ua}<Ja} • Thus, Oq 



commutes with the local POVMs Eq. (32) on all v £ C, 



Second, consider the neighboring domains C p 



Or 



by construction. Oc thus commutes 



with the local POVMs F v , a for all v £ C M and for all fj,. 

Therefore, Oc a is in the stabilizer of |^(^4)). Therefore, 
Oc a is an encoded operator w.r. t. the code in Table [TTJ 
and we need to figure out which one. (1) Central ver- 
tex C c : Oc c \c c is an encoded operator on G c , Oc c \c c G 
{±7, ±X, ±F, ±Z}. S_ince_0 Cc |, ? a$ for any l£ C c ,hy 
Table|nl Cc \ Cc ^±I,±Z. Thus, O c Jc c £ {±X,±Y}. 
(2) Neighboring vertices C M : By Table [n] a^l = ±Z, for 
any I £ Cn. Thus, 0, 



= ±Z 



1^1 



Now observe that 



Z = /, and that, this justifies the above prescription in 
constructing the graph G{A). Using the adjacency ma- 
trix Ag(a)i we have mod 2 = [Aq^] c ^ and hence 



O, 



c c \c u 



_j_T^G(A)] 



B a = {(| + ++)±|-— ))A/2}, 
6, = {(|i,i,i>± I -*>-*, -*»/V2} 
6 2 = {(|000)±|lll))A/2}. 



(35) 



These measurements map the symmetric subspace of the 
three-qubit states into itself and they can therefore be 
performed on the physical spin 3/2 systems. 

Denote by Sc and Sc\ v the code stabilizer on the do- 
main C £ V(G(A)) and on the reduced domain C\v, re- 
spectively. Using standard stabilizer techniques [38| it 
can be shown that the measurement Eq. ( 35 ) has the 
following effect on the encoding 



±X 



C\v, 



J C\v 



(36) 



The measurement ( 35 1 thus removes from C by one lattice 
site v £ V(C). We repeat the procedure until only one 
site, w, remains in C, for each C £ V(G(A)). In this 



way, S c 

MA)) - 



Ulo, 



X 



c 



±x 



{w}> 



Z 



Thus, 



Z \G(A)) =: \G(A)), where U\ oc is a local 
unitary, and the encoding in Table [TT] has now shrunk to 
one site of C per encoded qubit, i.e. to three auxiliary 
qubits. 

To complete the computation, the remaining encoded 
qubits are measured individually. Again, the measure- 
ment of an encoded qubit on a site w £ C is an operation 
on the symmetric subspace of three auxiliary qubits at 
w, and can thus be realized as a measurement on the 
equivalent physical spin 3/2. □ 



RANDOM GRAPH STATES, PERCOLATION 
AND 2D CLUSTER STATES 



Whether or not typical graph states \G(A)) are univer- 
sal resources hinges solely on the connectivity properties 
of G(A), and is thus a percolation problem [JD]. We test 
whether, for typical graphs G(A), 
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CI The size of the largest domain scales at most loga- 
rithmically with the total number of sites |V(£)|. 

C2 Let S C C be a rectangle of size I x 21 (21 x I). Then, 
a path through G\s from the left to the right (top 
to bottom) exists with probability approaching 1 in 
the limit of large C. 

Note that Condition CI is obeyed whenever the domains 
are microscopic, i.e., their size distribution is independent 
of \C\ in the limit of large C. Then, the size of the largest 
domain scales logarithmically in |£| [40 . Condition C2 
ensures that the system is in the percolating phase. 

Together with planarity, which holds for all graphs 
G(A) by construction, the conditions CI and C2 are suf- 
ficient for the reduction of the random graph state to 
a standard universal cluster state. The proof given be- 
low extends a similar result already established for site 
percolation on a square lattice [H]. The physical intu- 
ition comes from percolation theory. In the percolating 
(or supercritical) phase, the spanning cluster contains a 
subgraph which is topologically equivalent to a coarse- 
grained two-dimensional lattice structure. This subgraph 
can be carved out and subsequently cleaned off all imper- 
fections by local Pauli measurements, leading to a perfect 
two-dimensional lattice. 



A. Reduction of \G(A)) to a 2D cluster state above 
the percolation threshold 

We define the distance distc (v, w) between two vertices 
v, w € V(C) as the minimum number of edges on a path 
between u and v, and consider two further properties of 
graphs G(A): 

CI' G(A) can be embedded in C such that the maxi- 
mum distance between the endpoints of an edge in 
E(G) scales at most logarithmically in |£|, 

C2' Let S C C be a rectangle of size 1x51 (5lxl). Then, 
a path through G\s from the left to the right (top 
to bottom) exists with probability approaching 1 in 
the limit of large C. 

Lemma 1 G(A) is planar for all POVM outcomes A. 
Property CI implies property CI' , and property C2 im- 
plies property C2' . 

Proof of Lemma^ Planarity: G(A) is obtained from 
the honeycomb lattice C, which is planar, by the graph 
rules Rl and R2. They only perform edge deletion and 
edge contraction, which preserve the planarity of C. G(A 
is thus planar for all A. 

Property CI': For any domain d C V(C), place the 
corresponding vertex v(d) £ V(G(A)) inside d such 
that the distance r(d) := max w ^d distc(w, v) is mini- 
mized. Then, r(d) < \d\/2. Now, consider two vertices 
v(di),v(d 2 ) £ V(G) connected by an edge e £ E(G). 
Then, the domains d\ , di are connected by a single edge 



in E(C). Thus, for any pair (v(dx), i>(g?2)) of vertices in 
G connected by an edge in E(G), distc (v(di),v(d 2 )) < 
(\di \ + |rf 2 |)/2 + 1. By property CI, this length scales at 
most logarithmically in \C\. 
Property C2': See Fig.[6ji. 

Lemma 2 Consider a planar graph G(Ac) embedded 
into the lattice £ of size A x A, satisfying the proper- 
ties CI' and C2' . Then, the graph state \G(Ac)) can 
be converted by local measurements to a two-dimensional 
cluster state of size A' x A', with A' ~ A/ log A. 

Both lemmata combined give the desired result: 

Theorem 2 Consider an AKLT state on a honeycomb 
lattice C converted into a random graph state \G(A)) 
by the POVM §3(fy. If for typical POVM outcomes A 
the corresponding graph G(A) satisfies the conditions CI 
and C2, then the AKLT state is a universal resource for 
MBQC. Furthermore, the computation requires at most 
a poly-logarithmic overhead compared to cluster states. 

Remark: The polylog bound to the overhead comes 
from bounding the average domain size by the maximum 
domain size, for technical reasons. The true overhead is 
expected to be constant. 

That the conditions CI and C2 are obeyed in the typ- 
ical case remains to be demonstrated. We show this nu- 
merically, as reported in the next section. 

Proof of Lemma [J| - main tool. We show that a graph 
state \G) satisfying the assumptions of Lemma [2] can be 
reduced to the cluster state on a two-dimensional square 
grid, by local Pauli measurement on a subset of its qubits. 
The 2D cluster state is already known to be universal 
[8]. Specifically, we use the following rules [39] for the 
manipulation of graph states 



(37a) 






(37b) 



(37c) 



Rule ( |3"7a| ): The effect of a o~ z -measurement at vertex 
a on the interaction graph is to remove a and all edges 
ending in a. 



Rule (37b |: The effect of a (7,,-measurement at vertex 



a on the interaction graph is to invert all edges in the 
neighborhood N(a) of a, and to remove a and all edges 
adjacent to a. 

Rule (37c): Consider three qubits on a line, where the 
middle qubit has exactly two neighbors. When the left 
and the middle qubit are measured in the o^-basis, the 
interaction graph changes as follows: The right vertex 
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a) 



b) 



FIG. 6. (a) Constructing a rectangle with a traversing path 
of aspect ratio 1:5, from such rectangles with aspect ratio 1:2, 
2:1. (b) Overlapping rectangles of size L x 5L, 5L x L. The 
union of their traversing paths yields the net V ■ (c) Pair of 
non-separated junctions. 



™ 1 V \ 


a ppSa 






















FIG. 7. Transforming jG(.4)} into a 2D cluster state, (a) 
Macroscopic view: Regions a, /?, A etc imposed on the graph 
G(A). (b) Graph with three-valent junctions, (c) Decorated 
2D grid, (d) 2D grid for the cluster state. 



inherits the neighbors of the left vertex. The left and the 
middle vertex plus all edges adjacent to them are deleted. 

Proof of Lemma [2] - Outline. We consider a graph G 
with properties CI', C2'. We impose a pattern of regions 
A, B, C, .. of rectangular shape and size I x 51 and 
51 x I, for sufficiently large I, on the plane into which G is 
embedded; See Fig. [6]d. Due to the percolation property, 
G has a net-shaped subgraph V, shown in Fig. [7^,. In the 
first step of the reduction all qubits in V(G)\V(P) are 
measured individually in the a z -basis. The graph thereby 
created is close to the one displayed in Fig. [Tfj. However, 
it may have ad dition al edges that cannot be removed by 
vertex deletion ( 37a ) alone. Su ch ed ges ar e rem oved by a 
combination of the graph rules ( 37a I and ( 37b I . Then, in 



two further steps, the graph state of Fig. [7b is converted 
to the 2D cluster state shown in Fig. [7|i. 

Step 1. - Conversion of the graph G to the graph of 
Fig.[7)D by local operations. First we show that the net V 
of paths shown in Fig. [7^, exists. Consider the overlapping 
rectangles A and B in Fig. ^p. By Property C2', A has 
a path V a running from top to bottom, and B a path 
Vb crossing from left to right. Since G is planar, Va and 
Vb must intersect in at least one vertex. The net V is 
defined to be the union of all such traversing paths (one 
per rectangle), with all ends removed that do not affect 
connectedness. V is shown in Fig. [7£i. 

Now, G is converted to G\v(v), by deleting all vertices 
from G which are not in V(V). Physically, the corre- 
sponding graph state |G|y(-p)) is obtained by measuring 
the qubits at all vertices v <G V(G)\V(V) in the eigenba- 
sis of a z , c.f. Eq. (37a). After that, ideally, all junctions 
of paths in G|y(-p) should be T-shaped, as in the graph 
of Fig. [TJd. However, in general they will not be. The 



quintessential (but not only) obstruction is 




Note that no further edges can be removed by vertex dele- 
tion (corresponding to ^-measurements), without dis- 
connecting the junction. Nonetheless, this obstruction is 
easily dealt with. The above ring junction is converted 
into a T-junction by a single measurement in the a y - 
basis, 




(38) 



However, we have to show that all possible obstructions 
to T-junctions can be removed. 

We begin with the wires, running from one junction 
Jl G V{G\ V (v)) t° another junction J R 6 V(G\ V (v))- 
They also have obstructions, for example 



J, 



L 



J, 



R 



The first goal is to remove all obstructions within each 
wire. We require that for a given wire, each belonging 
vertex, except and Jr, has a unique neighbor to the 
right and a unique neighbor to the left in the wire, to 
which it is connected by an edge. Jl (Jr.) only has a 
unique right but no left (a unique left but no right) neigh- 
bor in the wire. The vertices are allowed to have edges 
with vertices in other wires. Those will be removed later. 



12 



Note that the vertices in the wire are left-right ordered, 
by order of appearance in the corresponding percolation 
path. Now, the obstructions are removed from any given 
wire by the following procedure. Starting with v = Jl, 
take a vertex v and identify its rightmost neighbor in the 
wire, w. Delete all vertices between v and w. Now set 
v := w, and repeat until w — Jr. 

In this way, each vertex in the wire remains with a 
single right neighbor in the wire. Therefore, the number 
of edges within the wire equals the number of vertices 
minus 1. Therefore, each vertex in the wire also has a 
unique left neighbor. 

At this stage we remain with the obstructions at junc- 
tions. First we show that we can treat the junctions 
individually, by choosing a sufficiently large length scale 
I for the size I x 51- (51 x 1-) rectangles. Two neighbor- 
ing junctions Jl, Jr have a distance of at least I, c.f. 
Fig. [BJd. They are separated if the configuration of edges 
shown in Fig. |6j: does never occur. It doesn't occur if 
I > 2|e| max , where |e| max is the maximum distance of an 
edge in E(G(A)). With CI' it thus suffices to choose 



logk. 



(39) 



Now we discuss an individual junction J. By construc- 
tion it joins three wires, Wl, Wc and Wr say. The 
obstructions are three sets of edges, Elr, Elc and Ecr- 
They connect vertices in Wl with vertices Wr, Wl 
and Wc, and Wc and Wr, respectively. By the choice 
Eq. (39) for I, the obstructions at different junctions are 



well separated from each other, and we can thus treat 
them individually. 

First, we remove the obstructions Elc an d Ecr, by 
the following procedure. We approach the junction at 
J on the wire Wc- Denote by v the first vertex in Wc 
which is the endpoint of an edge e £ Elc U Ecr- By 
rule (37a), delete all vertices in Wc between v and J, 
excluding v and J. Then, there arise three cases. (1) 
v is connected to a single vertex in Wl U Wr. (2) v is 
connected to exactly two vertices w\,W2 & Wl U Wr, and 
wi, W2 are neighbors in Wl U Wr. (3) v is connected to 
exactly two vertices wi,W2 € Wl U Wr, and 101,102 are 
not neighbors in WlUWr; or v is connected to more than 
two vertices in Wl U Wr. Graphically, the cases look as 
follows (the obstructing edges Elr are not relevant in 
the present sub-step, and are not shown), 



(3) 



Case (1): The vertex w is taken as the new junction cen- 
ter J, and the edge (v, w) is included into Wc- Thereby, 
the obstructing edge sets Elc and Ecr are removed. 
Case (2): The qubit on vertex v is measured in the a v - 
basis; c.f. Eq. (38). Thereby, case (2) is reduced to case 
(1). Case (3): Denote the leftmost (rightmost) neighbor 




of v in Wl U Wr by wl (wr). All vertices in Wl U Wr 
between wl and wr are deleted, by <7 2 -measurement on 
the corresponding qubits. Thereby, case (3) is reduced 
to case (1). 

In the above procedure, the center J of the junction 
may have shifted within Wl U Wr. Consequently, Wl 
(to the left of J), Wr (to the right of J) and Elr are 
modified. Due to the shift in the location of J, edges that 
were in Elr may have become obstructing edges inter- 
nal to Wl or to Wr. They are removed by re-running 
the previous procedure for removing obstructing edges 
internal to the wires. There are no new edges in Elr- 

By the above procedure, we have created a junction 
of three wires Wl, Wc and Wr which are free of all 
obstructions except for the set Elr- These edges are 
now removed by approaching the new junction center J 
from the wire Wl and repeating the previous procedure. 

Step 2 - Creating the decorated lattice graph of Fig. [7b. 
Consider a ring-shaped segment of the graph in Fig. [7p, 
with the four belonging T-junctions a, j3, 7 and 5. The 
qubits on all vertices on the path between a and /3 are 
measured in the cr z -basis. The corresponding vertices 
are thereby removed, c.f. Eq. (37a). Regarding the path 
between a and 5, we require that a and j5 are not neigh- 
bors. This can always be arranged by starting with a 
sufficiently large scale I for the I x 5/-rectangles. If the 
the number of vertices that lie on the path between a and 
6 is even (but > 0), then the qubit on the vertex next to 
a is measured in the <7j,-basis. In this way, the number 



of vertices between a and 5 becomes odd, c.f. Eq. (37b). 
Now, a and the vertex next to it are measured in the 
cr^-basis. By Eq. (37c), a moves two vertices closer to S 
This procedure is repeated until a and /3 are merged into 
a single vertex a' . Then, in an analogous manner, a' is 
merged with 7, and with /3. Thereby, the ring of four 
T-junctions is converted into a single vertex of degree 4. 

Step 3 - Creating the square lattice graph of Fig. [TJl. 
By the same method as in Step 2, the line segments 
between the vertices of degree 4 are contracted. This 
creates a two-dimensional cluster state on a square lat- 
tice, which is known to be a universal resource for 
measurement-based quantum computation [8]. Since 
only local measurements were used in the reduction, the 
original graph state \G) is universal as well. 

Overhead: The bottleneck of the construction is to 
guarantee that all junctions can be treated individually, 
which requires L ~ log A, c.f. Eq. (39). As displayed in 
Fig. [6b-, a qubit in the created planar cluster state claims 
an area of size 81 x 81 on the original honeycomb lattice 
C, and thus A' ~ A/ log A, as claimed. □ 

Our proof thus generalizes the results in Ref. [41] to 
random graphs. This extends the set of cluster- type uni- 
versal states to more general 2D random graph states and 
beyond regular lattices [TT] , 
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VI. MONTE CARLO SIMULATIONS 

We give the recipe for performing Monte Carlo 
simulations and present some results. 

(1) First, we randomly assign every site on the hon- 
eycomb lattice to be either x, y or z-type with equal 
probability. 

(2) Second, we use the Metropolis method to sam- 
ple typical configurations. For each site we attempt 
to flip the type to one of the other two with equal 
probability. Accept the flip with a probability 
Accept = min{l,2^'H£'Hi'l+|£|} ) wncrc \V\ and 

\£\ denote the number of domains and inter domain 
edges (before the modulo-2 operation on inter-domain 
edges; see Fig. 2 of main text), respectively before 
the flip, and similarly \V'\ and \£'\ for the flipped 
configuration. The counting of \V\ and \£\, etc. is done 
via a generalized Hoshen-Kopelman algorithm [42]. For 
the proof of the probability ratio, see Sec. |VI A| 

(3) After many flipping events, we measure the prop- 
erties regarding the graph structure for the domains 
and study their percolation properties upon deleting 
edges. For the percolation, we cut open the lattice and 
investigate the percolation threshold for the typical 
random graphs from the Metropolis sampling. 



A. Evaluation of probability ratio 

In this section, we shall explain the transition proba- 
bility ratio: Accept = min|l,2l v 'l-l £ 'l-l v l+l fi l}, which 

arises from the probability for POVM outcomes {a v } 
being p({a v } ~ 2l y l~l £ l, where |V| and \£\ denote the 
number of domains and inter-domain edges (before the 
modulo-2 operation). The proof is very similar to that 
in ID. For convenience, we shall use spin-3/2 representa- 
tion of the AKLT state. The local mapping from three 
virtual qubits to one spin-3/2 is 

Pv = |1)(000| + |2)(111| + |3)<W| + \4)(W\, (40) 

where we have simplified the notation for the spin-3/2 
basis states: |1) = |3/2,3/2), |2) = |3/2, —3/2), |3) = 
|3/2, l/2)_and |4) = |3/2, —1/2). Moreover, |000), | 111) , 
\W) and \ W) constitute the basis states for the symmetric 
subspace of three spin-1/2 particles. The AKLT state can 
then be expressed as 



/AKLT 



®a n 



(41) 



where \<j>) e is the singlet state ( 1 01) — |10)) Uii ^. for the 
edge e = (u, v) and i,j specify the virtual qubit in the 
respective vertex. 

The POVM that reduces the spin-3/2 AKLT to a spin- 
1 /2 graph state consists of elements = F^F^ such that 



1 = E x + E y + E z , with 
Fz = Ft = 
Fl = 



F„ 



Ft 
y 



y|(|l)(l| + |2)(2|) = -^(5 2 2 -i), (42) 

(43) 
), (44) 



= (|a)<o| + |6)<6|) = ^(^ 



l 



i 



3 (|a)H + |/3)(/3|) = -^(5 y 2 



where we have also expressed F"s in terms of the cor- 
responding spin operators. The other four states other 
than |1) and |2) are 

\a) = \S X = 3/2) = i(|l) + |2) + V3|3) + >/3|4)) (45) 

\b) = \S X =-3/2) = J=(|l) - |2) - s/3\3) + V3\4)) (46) 

\a) = \S y = 3/2) = -L(|l) - i\2) + zV3\S) - >/3|4)) (47) 

|/3) = \S V =-3/2) = -L(|l) + i\2) - zV3|3) - s/3\4)).(48) 

They correspond to the four virtual three-spin-1/2 states 

(in addition to |000) and | 111)) | + ++), | ), \ i i i) 

and | — i). 

While the outcome of POVM constructed above at 
each site is random (either x, y or z), outcomes at dif- 
ferent sites may be correlated. For a particular set of 
outcomes {a v } at sites {?;}, the resultant state is trans- 
formed to the following un-normalized state 



W) 



AKLT, 



with the probability being 

P{a„} = W)/(V#)AKLT. 



(49) 



(50) 



As F 's are proportional to projectors, in evaluating the 
relative probability for two sets of outcome {a v } and 
{b v }, one has 

P{a v }/P{b v } = (^l(^)^,aJV')AKLT/(V'l(^)^.fcJV ; )AKLT, 

V V 

(51) 

where we have used F^ a ~ F v a . In order to evaluate the 
probability ratio for two different sets of configuration, 
we first note that 

F,p~|o)<+ + +| + |6)< 1 (52) 

FyP~\a){iii\ + \P){-i-i-i\ (53) 
F Z P ~ 1 1) (000 1 + |2)(111|. (54) 

The spin-3/2 state is transformed by F ViClv to an ef- 
fective spin-1/2 one, with the two levels being labeled by 
(a, b), (a, /3), or (1,2), depending on which F is applied. 
The probability ps av \ is essentially obtained by summing 
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the square norm of the coefficients for all possible spin- 
1/2 constituent basis states (e.g. \ab + Oi . . . ) is a basis 
state). First we need to know how many different con- 
stituent states, and the number is related to how many ef- 
fective spin-1/2 particles we have. For the sites that have 
same type of outcome (x, y or z), they basically form a 
superposition of two Neel-like states, thereby correspond- 
ing to an effective spin-1/2 particle. This can be seen 
from the valence-bond picture that, e.g., for r,s € {0, 1} 
we have (rs|01 — 10) = ±t5 ri i_ s . On the other hand, 
for r e {0,1} and s € {+,-}, (rs|01 - 10) = ±1/V2, 
which is 1/ v2 smaller than if r and s are the indices 
in the same basis. This means that all four combina- 
tions {0+, 0—, 1+, 1 — } occur with equal amplitude up to 
a phase. (Similar consideration applies to other combina- 
tions of bases.) Therefore, the number of effective spin- 
1/2 particles is given by the number of domains, which 
we label by \V\. Notice that we have assumed that any 
domain does not contain a cycle with odd number of orig- 
inal sites, as no Neel state can be supported on such a 
cycle (or loop). Configurations with domains that con- 
tain a cycle with odd number need to be removed. For- 
tunately, as the honeycomb lattice is bi-partite, therefore 
any cycle must contain even number of sites and we do 
not to deal with the above complication. 

What about the amplitude for each spin configuration? 
Furthermore, what is the probability of obtaining a par- 
ticular set of outcome {a„}? We have seen that for each 
inter-domain edge there is a contribution to a factor of 
1 / y2 in the amplitude (as the end sites of the edge corre- 
spond to different types). Thus, the amplitude for each 
spin configuration gives an overall value (omitting the 
phase factor) of 2~l £ l/ 2 and hence a probability weight 
2~' £ ', where \£\ counts the number of inter-domain edges. 
As there are 2'^' such configurations, we have the norm 
square of the resultant spin-1/2 state being proportional 
to p ~ 2l^l _ l £ l . For convenience, we have assume the lat- 
tice is periodic, but the argument holds for open bound- 
ary condition in which the spin-3/2's at the boundary 
are either (1) suitably linked by to one another, preserv- 
ing the trivalence or (2) terminated by spin-l/2's. In the 
Appendices, we have provided an alternative derivation 
of the probability expression. 



B. Discussions of simulation results 

We have analyzed lattices of size up to 200 x 200 
sites. As shown in Fig. [HJ the size dependence of av- 
erage vertex number, average edge number, and aver- 
age Betti number B (43) of the random graphs formed 
by domains relative to the original lattice size behaves 
as follows: \V\ = 0.495(2)L 2 , \E\ = 0.872(4)L 2 , and 
B = 0.377(2)L 2 , where L is related to the total num- 
ber of sites in the original honeycomb lattice N = L x L. 
This shows that the typical random graph of the graph 
state retains macroscopic number of vertices, edges, and 
cycles, giving strong evidence that the state is a universal 
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FIG. 8. (color online) Average vertex (or domain) number, 
average edge number, and average Betti number in the typical 
random graphs original lattice site vs. L. The total number 
of sites is N = L 2 . This shows the number of domains, the 
number of interdomain Ising interaction, and the number of 
independent loops in the resultant graph all scale with the 
system size of the original honeycomb lattice. 



resource. Figure [9] shows the average degree of a vertex 
vs. inverse system length 1/L for the random graphs, 
as well as the average numbers of the original sites con- 
tained in a typical domain. The average vertex degree 
extrapolates to d w 3.52(1) for the infinite system. This 
compares to 4 for the square lattice and 3 for the honey- 
comb lattice. 

In order to show the stability of the random graph, we 
investigate how robust it is upon, e.g., deleting vertices 
(or edges) probabilistically, i.e., performing the site (or 
bond) percolation simulations. As shown in Fig. [TO^ i, it 
requires the probability of deleting vertices to be as high 
Pdeiete = 0.33(1) (i.e., percolation threshold p c = 0.67(1)) 
in order to destroy the spanning property of the graph. 
This lies between the site percolation thresholds « 0.592 
of the square lattice and ~ 0.697 of the honeycomb lat- 
tice. For bond percolation as shown in Fig.|10|b, it takes a 
probability of Pdciotc = 0.43(1) (i.e., percolation thresh- 
old p c = 0.57(1)) to destroy the spanning property of 
the graph. Again, this threshold lies between that of 
the square lattice (1/2) and that of the honeycomb lat- 
tice (« 0.652). This shows that there exists many paths 
(proportional to the system's linear size) on the random 
graphs that can be used to simulate one-qubit unitary 
gates on as many logical qubits and entangling opera- 
tions among them. We remark that percolation argu- 
ment was previously employed by Kicling, Rudolph, and 
Eisert in establishing the universality of using nondeter- 
ministic gates to construct a universal cluster state [46] . 

Let us also examine the two conditions listed in 
Sec.lVZI 
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FIG. 9. (color online) Average domain size (i.e., number of 
original sites in a domain), average width of domain size dis- 
tribution, and average degree of a vertex in the typical random 
graphs vs. L, where L 2 is the total number of sites in the hon- 
eycomb lattice. For better discernibility of the two lower sets 
of data, we suppress the errorbars for one of them. This set 
of data was first shown in Ref. [34], and we have reproduced 
it here for the sake of completeness. 



Condition CI. For all POVM outcomes sampled from, 
the size of the largest domain was never macroscopic and 
it can at best be logarithmic in the original system size; 
see Fig. 11 The average number of sites v £ V(C) con- 
tained in a typical domain, when extrapolated to the 



(b) 



infinite system, is 2.02(1); see Fig. [9| Our numerical sim- 
ulations thus show that condition CI holds. 

Condition C2. For all of the POVM outcomes sampled, 
a horizontal and a vertical traversing path through the 
resulting graphs G(A) always existed (without deleting 
any vertex or edge). Our numerical simulations show 
that our random graph arc deep in the supercritical phase 
and thus condition C2 holds. 

In addition, a necessary condition for the computa- 
tional universality of the graph states \G(A)} is that typ- 
ical graphs G(A) are not close to trees, because MBQC 
on tree-like graphs can be efficiently classically simulated 
|44j . For typical graphs G(A), we find that the Betti 
number (which is zero for trees) is proportional to the 
size of the initial honeycomb graph, with B = 0.377(2)jV. 

Robustness. We now quantify how deep typical graphs 
G(A) are in the connected phase of the percolation tran- 
sition. A first measure is the average vertex degree. A 
heuristic argument based on a branching process suggests 
that a graph has a macroscopic connected component 
whenever the average vertex degree is d > 2. This crite- 
rion is exact for random graphs of uniform degree |40j . It 
also holds surprisingly well for lattice graphs [45], which 
are the least random. In our case, the typical graphs 
G{A) have an average degree of 3.52, suggesting that the 
system is deep in the connected phase, which is confirmed 
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FIG. 10. (color online) Percolation study of the graph formed 
by the domains: probability of a spanning cluster p c i ust cr vs. 
the probability to delete a vertex (top panel a) or an edge 
(bottom panel b) pdeictc- The threshold for destroying the 
spanning cluster is around Pdciotc ~ 0.33 in deleting vertices 
and Pdcictc ~ 0.43 in deleting edges. This shows that the 
graph without deleting any vertex or edge is deep in the per- 
colated (i.e., connected) phase. We note that the result of site 
percolation is reproduced from Ref. [33], and the additional 
bond percolation result presented here is consistent with the 
picture that the typical random graphs lie somewhat between 
the honeycomb and the square lattice. 



by the percolation simulations. The existence of finite 
percolation thresholds (for both site and bond percola- 
tion) discussed earlier further supports the robustness of 
the connectedness. 



16 



jections: 

G^| 7 f)( 7 f | + |7f 4 )<7® 4 |. (56) 
It can be checked that 

4 

YA £ ^- + ^E Gt * G *=P5, (57) 
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FIG. 11. (color online) The largest domain size in the typical 
graphs vs. L, with N — L 2 being the total number of sites. 
The fitted curve to the largest domain size is 3.337 ln(AT) — 
5.566. The result is reproduced from Ref. [34] . 



where Ps is the projection operator onto the symmetric 
subspace of four qubits [23]. However, such a generalized 
measurement would yield four additional pairs of states 
{ik, 7fc} which are not mutually unbiased to one another 
nor to the eigenstates of the three Pauli operators. Due 
to this complication, whether the 2D AKLT state on the 
square lattice is universal for MQBC remains open. 
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VII. CONCLUDING REMARKS 

We investigated the measurement-based quantum 
computation on the AKLT states. First we provided 
an alternative proof that the ID spin-1 AKLT state can 
be used to simulate arbitrary one-qubit unitary gates. 
We extended the same formalism and demonstrated that 
the spin-3/2 AKLT state on a two-dimensional honey- 
comb lattice is a universal resource for measurement- 
based quantum computation by showing that a 2D clus- 
ter state can be distilled by local operations. Along the 
way, we connected the quantum computational universal- 
ity of 2D random graph states to their percolation prop- 
erty and showed those 2D graph states whose graphs are 
in the supercritical phase are indeed universal resources 
for MBQC. 

The key ingredient that has enabled our proof of com- 
putational universality for the (spin-3/2) 2D AKLT state 
on the honeycomb lattice is the generalized measurement 
in Eq. Q. How about the case of (spin-2) 2D AKLT 
state on the square lattice or any other lattices beyond 
trivalence is universal for MQBC. The AKLT spin-2 par- 
ticle can be regarded as four virtual qubits in the sym- 
metric subspace. Hence a naive extension of the POVM 
prompts us to consider the follow operators: 



F z = (|0 )(0 | + 11 }(1 |) 



F„ 



H ){- 



_®4w_ t 



(55a) 
(55b) 
(55c) 



Unfortunately, ^2 a is not proportional to the pro- 

jection onto the symmetric subspace. However, we can 
consider additionally the four states | 7 fc) (k = 1,..,4) 
such that their Bloch vectors point in the four diago- 
nal directions of a cube, i.e., (1, 1, 1)/V3, (— 1,1, 1)/v3j 
(— 1, — 1, 1)/V3, and (1, — 1, l)/v3, respectively. To- 
gether with their corresponding conjugate states \~/f.) 
having opposite vectors, we have four other sets of pro- 



Appendix A: Calculation of probability of a 
particular POVM outcome using 
Arovas-Auerbach-Haldane techniques 

In this appendix we provide an alternative formulation 
to the calculation of POVM outcome probability. This 
formalism has the potential of being applicable to a more 
general case. We give only the important ingredients 
here. 

Arovas, Auerbach and Haldane (AAH) [57J show how 
to represent arbitrary AKLT states as Boltzmann weights 
for nearest neighbour statistical mechanical models in 
the same spatial dimension as the quantum problem and 
how to represent calculations of equal time ground state 
expectation values classically. We are interested in two 
cases: the 5 = 1 one-dimensional case and the 5 = 3/2 
honeycomb lattice case. 

In both cases the operators of interest are proportional 
to projection operators onto maximal \S Z \: 



F„ 



2, (5 = 1) 
l/4]/%/6, (5 



3/2), 



(Al) 
(A2) 



where v = x, y or z and, for convenience, we have rescaled 
the prefactor in the definition of -F"s. For general spin 
S the operators S a are represented first in terms of 
Schwinger bosons, a, a' ', b, w, then in terms of co- 
ordinates and derivatives u, v, d u , d v acting on homo- 
geneous polynomials of 0(25). The operator (5 2 ) 2 is: 

(S z ) 2 = (l/4)(ot - tfbf = (l/4)(flU - d v v) 2 

= (1/4) (c\V + d v v 2 - 2d u d v uv - d u u - d v v). (A3) 

We now use the prescription of Arovas, Auerbach and 
Haldane: 

(^|a*4«*+'V-''|v> 



k+l+l 



n ( 2s 



'» (A4) 
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for any states and \ip) in the spin 5 Hilbert space. 



To prove Eq. ( A4 ) , note that a complete set of states 
for the spin 5 Hilbert space is given by u s + m v s - m which 
are eigenstates of 5 2 with eigenvalue m = — 5, — 5 + 
1, . . . 5. To prove Eq. (A4| for j — 0, we wish to prove: 



d 2 Q,u* s+m v* s - m dtdlu k+s+m v l+s - m 



k+l+l 1 

J] (25 + r) d 2 n\u\^ s+m+ ^\v\ 2 ^- m+l K(A5) 

. r=2 J "* 

To prove this, we use the identity: 



i M = / d 2 n\ u \ 2p \v\ 2 i 



2ir{l/2) p+q J dx(l + x) p {l-x) g 
4irp\q\/(p + q + iy. 



(A6) 



Thus the left hand side of Eq. ( A5 ) may be written 



TUC _ (S + m + k)\(S-m + iy. T 

LHS ~ (S + m)l (S-m)\ Is + m ' s ~ m 



= An 



(S + m + k)l(S-m + l)\ 



(25 + 1)! 



k+l+l 



n ( 25+r ) 



Is 



+m+k,S-m+l 



(A7) 



which is the RHS. Furthermore, all off-diagonal matrix 
elements vanish for both the left and right hand side of 
the identity in Eq. (A4) for j = 0. That follows since 



rp^tpmi oc e l ( m ~"W where <j> is the azimuthal angle for 
the integration over the sphere. Neither inserting the op- 
erator on the left hand side of Eq. ( A4 ) nor multiplying 



by the function on the right hand side changes this az- 
imuthal angle dependence, implying vanishing integrals. 
While it may appear that this proof only holds for j = 



in Eq. (A4) it actually covers the case of general j. In 
general, Eq. (A5| gives the (ip m \ . ■ ■ |^ ro _j) matrix ele- 



ments of the identity in Eq. ( A4 ) , which are the only 



non-zero matrix elements. Furthermore the identity im- 
mediately generalizes to an arbitrary product on different 
lattice sites: 



ki+h + l 

I] n (25 + m) 

i m—2 



-mil 



*L hi-\-ji L — ii 
U 4 % Va % Ui J t>- 



(A8) 



since we may simply extend the above argument to the 
basis states n» u i mi vf~ mi for which the matrix ele- 
ments simply factorize. Since we have proved this iden- 
tity for a complete set of states it follows for any states 
\ip') in the spin-5 Hilbert space, including the AKLT 
states. 



Eq. (A4) gives: 

{S z ) 2 = (l/4)(25 + 2)(25 + 3)[M 4 + \v\ 
-(1/4)(2S + 2)(M 2 + H 2 ) 
= (1/4)(2S + 2)(2S + 3)[M 2 - M 

-(i/2)(s + i)(M 2 + M 2 ). 



2M 2 M : 



2l2 



(A9) 



Using u = cos(0/2)e^ /2 , v = sin(6'/2)e-^ /2 , where 6 and 
(j) are the polar and azimuthal angle on the unit sphere, 
this becomes: 



(5 



2\2 



(l/4)(25+2)(25+3)(^) 2 -(l/2)(5+l) (A10) 



where Q z — cos 9 is the projection of the unit vector onto 
the z-axis. A simple explicit calculation similar to this 
one shows that, for v = x, y or z: 

(5„) 2 = (l/4)(2S+2)(2S+3)(rn 2 -(l/2)(S+l) (All) 

as expected by SO(3) symmetry. This is a somewhat 
surprising formula in that the classical quantities are not 
positive semi-definite. Note that this formula is valid 
independent of the wave-function. The projection oper- 
ators thus become: 

F v = [5(fi") 2 — 1]/V2, (5 = 1) (A12) 
= y/3/8[5((V) 2 - 1], (5 = 3/2). (A13) 

Remarkably the projection operators are the same for 
5 = 1 and 3/2 up to an unimportant normalization fac- 
tor. 

The AKLT state can be written, in Schwinger boson 
notation as: 



) aklt = II( a H ~~ a] &t) | vacuum) 



corresponding to 

ip(v>i,Vi 



n ( Uiv i 

(id) 



(A14) 



(A15) 



where the product is over all pairs of neighboring sites 
. The square of the wave-function is: 



\i>(Ui,Vi)\ 2 cx 

(id) 



(A16) 



Actually, we need to be more precise about boundary 
conditions here. These details will be discussed below. 
Using the form of the AKLT state we wish to calculate: 



dfli[5(Q^) 2 - 1] Yl [1 - • Oj] 

(J,k) 

(A17) 

where Z is the same integral without the [5(f2"*) 2 — 1] fac- 
tors and Af s = (1/2)" for 5 = 1 and (1/3)" for 5 = 3/2. 
Note that the inserted operators are F u /y2 for the 5=1 
case and iJ2/3F u for the 5 = 3/2 case, normalized so 
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(a) (b) 
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FIG. 12. Illustration of loops, (a) A domain consisting 
of 8 vertices, 8 edges and 1 face. The numbers indicate one 
possible order in which the domain could be grown. When 
the first vertex is added, V = 1, E — F — 0. As vertices 2 
to 7 are added both V and E increase by 1. When the 8^ 
vertex is added, V and F increase by 1 and E increases by 2. 
(b) In this domain, which contains F = 3 faces, two of them 
are selected, as indicated by shading. The corresponding set 
of loops is a single loop surrounding the two adjacent faces 
as indicated by heavy lines. This domain contains V = 16 
vertices and E = 18 edges, obeying F = E — V + 1. (c) This 
is the same domain as in the previous figure (b) but a different 
subset of faces is selected. Now the set of loops consists of 2 
loops as indicated by heavy lines. 



that the sum over v gives the identity operator, ensuring 
the proper normalization of the probability distribution. 
In both cases we can evaluate this by multiplying out 

Carrying out this we arrive at the same conclusion of 
the probability expressions for ID chain and 2D honey- 
comb cases as before, as we show below. 



Appendix B: S=l, 1 dimension 

We first consider the ID S=l case as a warm-up. 



1. Open Boundary Conditions 

Consider a chain of n spin-l's on sites i — 1, 2, ... n 
with 2 additional S=l/2's at sites and n+1 to remove 
the "dangling bonds" . Then the AKLT ground state is: 

n 

\^)^\{{a\b\ +1 -a\ +1 b^\0). (Bl) 

i=0 



Thus: 

n 

|^o| a = n[i-n ( -A i+1 ]. (B2) 

We only make projective measurements on the sites con- 
taining spin-l's, at 1, 2, . . . n. In this case we may replace 
each factor 1 — f2j • f2j by 1 because all other terms in the 
expansion contain a single power of one or more vector 
and thus give zero after integrating over Oj. Using: 

((fi°) 2 ) = (l/3)((tt) 2 ) = 1/3, ((50") 2 -l) = 2/3, (B3) 

we obtain a constant: 

pai a 2 ...a n = (j^n ( B4 ) 

independent of the Oj's. 

2. Periodic Boundary Conditions 

Now we consider n sites, all with spin-l's and cou- 
ple site n to site 1. This is a useful warm up for the 
2D case because there is now one closed loop, i.e., one 
other term in the expansion can give a non-zero integral: 
( — 1)™ n"=i ' Now we need the integral: 

J dn[5{n a f - i]n b n c . (B5) 

Clearly this vanishes unless b — c. Note that: 

((IF) 4 ) = (cos 4 fl) = (1/2) ^ dxx 4 = 1/5 (B6) 

and 

<(n*) 2 (lF) 2 ) = fj g cos 2 0sin 2 f(l - xr 

= 1/15. (B7) 
Thus we obtain the remarkable identity: 

([5(fT) 2 - l]il b n c ) = -5 bc 5 ab . (B8) 
3 

The integral vanishes unless a = b = c in which case 
it has the same value as ((5f2 a ) 2 — 1). The product 
( — 1)™ Y[7=i &i ' f^t+i contains sums over n indices. How- 
ever, for the integrals to be non-zero all cti indices must 
equal each other. In this case the multiple integral has ex- 
actly the same value as when the product is not present, 
giving: 

pa 1 a 2 ...a n _ /i /o\n [j; + (~ Q2 ^0203 • ■ • 

- I. I*) l + (-l)'»(l/3)'»- 1 

Here we have used the fact that the partition function 
also obtains a contribution from (—1)™ n"=i ^i'^i+i 
ing the second term in the denominator and the reason 
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that it is three times larger than 1 /3" is because of three 
possibilities a = x,y, z, or more precisely, 



(a) 



/n 
Y[d 2 n t n t -n l+1 = (i/3)"-\ 
i=i 



(BIO) 



where f2 n +i = S~2i . Thus, p°i a 3— a n j g near ]y constant 
again except that in the one case where a\ = a-i = . . . a n 
it is twice as big if n is even or zero if n is odd. This result 
agrees precisely with that obtained by other methods in 
sub-section Sec. Ill Dl 



Appendix C: S=3/2, 2 dimensions 

1. Open Boundary Conditions 

Consider an arbitrary finite segment of a honeycomb 
lattice, consisting of n spins; it could have zig-zag and 
armchair edges or disordered ones, for example. Spins 
on the boundary will generally be coupled to cither 2 or 
3 other spins - 2 for a zig-zag edge and 3 for an armchair 
edge, for example. In all cases where a boundary spin 
is only coupled to 2 other spins, couple it to a boundary 
S=l/2 spin. Let the total number of spins, including the 
S=l/2 spins on the boundary be M. Then the square of 
the AKLT ground state is: 



(b) 



\1>o\ 2 =1l[l-tk-&j]. 



(CI) 



The product is over all nearest neighbours, as usual in- 
cluding both S=3/2 and S=l/2 spins. We only do the 
POVM on the S=3/2 spins. Since each of the boundary 
S=l/2 spins couples to only one other (S=3/2) spin, we 
may replace [1 — fij • Clj] by 1 for each factor involving an 
S=l/2 spin in calculating p a i a 2 - a n ^ Following the above 
reasoning, when we take the 1 term in the expansion of 



1 



(4tt) m (1/3)" (2/3) 



(C2) 



There will be many additional terms in this case, unlike 
the D=l case. Each additional term must correspond 
to a set of closed loops on the lattice, with zero or two 
lines entering each of the S=3/2 sites. These loops never 
involve the S=l/2 boundary sites. These loops can never 
cross each other but we can have loops inside loops. Such 
a contribution only exists when all the a^s for sites on a 
given loop have the same value. Each such term makes 
an equal contribution to p o i°2-" «.. Thus we simply need 
to calculate the number of sets of closed loops with equal 
dj's for a given configuration ai, ■ ■ ■ a n - 

To do this it is convenient to divide up all sites on the 
lattice into domains such that ai has the same value for 
all sites in a domain and all sites in a domain are the 
nearest neighbor of at least one other site in the domain. 



(c) 



FIG. 13. Illustration of loops in the case of the peri- 
odic boundary condition. A honeycomb lattice with periodic 
boundary conditions, drawn as a "brick wall" lattice for con- 
venience, (a) A topological loop is shown. The arrows indi- 
cate an edge in the domain between vertices at the left and 
right hand sides of the lattice, (b) The set of loops (one loop 
in this case) is shown for a domain containing one topological 
loop and one face sharing on edge with the topological loop, 
corresponding to the case in which the face and the topolog- 
ical loop are chosen. (c)A domain containing two topological 
loops, with W = 2. It can be seen that there are 3 possible 
loops, corresponding to 4 sets of loops, and zero faces. 



(Here sites refers to the sites with S=3/2 spins only.) The 
number of sites in a domain can range from 1 to n, in 
principle, although we expect that typical domains are 
microscopic. We draw a line between all nearest neigh- 
bors in each domain. We may identify a unique number 
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of faces with each domain, Fi and a total number of faces, 
F = Fi with a given configuration. A face, inside a 
domain, is an elementary hexagon which is completely 
surrounded by 6 edges belonging to that domain. Thus, 
it is impossible to move from the interior of a face to its 
exterior (either inside the domain or not) without cross- 
ing an edge belonging to its domain. The total number 
of sets of closed loops, Nl is then simply 



N L = 2' 



(C3) 



This follows because there is a unique set of loops which 
surrounds any subset of the faces. See Fig. [12] 

The i^ 1 domain will also have a number of edges, Ei 
and a number of vertices, Vi. It can be seen that: 



Fi = Ei -Vi + 1. 



(C4) 



This can be seen by induction, growing the domain site by 
site, always adding new sites which are nearest neighbors 
of at least one previous site. After drawing the first site, 
Vi = 1 and Ei = F t = 0, so Eq. fCl} is obeyed. When 



we add the next site, we increase both E{ and Vi by 1, 
without changing Fi. This goes on for a while but we 
may eventually add a site which is the nearest neighbor 
of 2 previous sites. At that step, V, increases by 1, E{ 
increases by 2 and Fi increases by 1 since we are then 
closing a loop, making a new face. Thus Eq. ( |C4[ ) remains 
true at each step as we grow the domain, completing the 
proof. 

Suppose we define a new, random lattice, by collaps- 
ing each domain down to one vertex, with an arbitrary 
number of edges, inherited from the original lattice, con- 
necting the various domains. Let V by the number of 
vertices of this random lattice and E' be the number 
of edges. (We ignore the S=l/2 boundary spins here.) 
Then: 



(C5) 



since Vi sites are reduced to 1 at the domain. Simi- 
larly if tie is the total number of edges connecting S=3/2 
spins in the original honeycomb lattice, then 



Thus: 



E' = n E -J2 E i 



(C6) 



pa ia2 ...a n ^ 2 V ( C7 ) 

the same result obtained in Sec. VI A| by another method. 



2. Periodic Boundary Conditions 

Now consider a honeycomb lattice of S=3/2's (no 
S=l/2's now) with periodic boundary conditions. This 
can be done in such a way that every spin has 3 nearest 
neighbors and we take the AKLT ground state. Similar 
to the D=l case, there can now be additional sets of loops 



because we can form loops that wrap around the torus 
but don't correspond to faces; see Fig. 13 If a domain 



wraps around the torus one way, but not the other, then 
the total number of sets of loops, corresponding to the 
domain is: 



Nr. 



vFi+l 



(C8) 



To see this choose an arbitrary "topological loop" within 
the domain going around the torus which doesn't encircle 
any faces. There are now 2 sets of loops corresponding 
to an arbitrary subset of faces, not using this topological 
loop and the arbitrary subsets of faces together with the 
topological loop. The construction of the set of loops cor- 
responding to the topological loop plus subset of faces is 
constructed by analogy with the above construction. In 
cases where none of the faces share edges with the topo- 
logical loop the set of loops corresponds to the topological 
loop plus the loops around the subset of faces. In cases 
where one or more faces shares an edge with the topolog- 
ical loop, the topological loop is modified to enclose each 
such face; see Fig. [T3p . 

Finally, it is possible to have 2 topological loops in 
a domain, going around the torus the two inequivalent 
ways; see Fig. 13 (In this case all other domains must 



be topologically trivial.) We can grow the domain ini- 
tially by drawing these 2 topological loops without any 
faces. At this stage there are 3 closed loops, going around 
the torus one way or the other or using all edges in the 
domain to go around both ways. Thus the number of 
sets of loops is 4 at this stage. After growing the entire 
domain we how have: 



N Li = 2 



Fi+2 



(C9) 



since each set of loops corresponds to a subset of faces, 
possibly combined with one of these three topological 
loops. In general, we may associate a winding number 
with each domain Wj = 0, 1 or 2 with: 



N Ll = 2 1 



(CIO) 



It can be seen that the number of vertices and edges 
in each domain obeys, in general: 



Fi + Wi=Ei-Vi + l. 



(Cll) 



This follows by induction, as we grow the domain. At 
the step when we complete the first topological loop we 
increase Ei by 2 but Vi by 1 and Fi by 0. Adding further 
faces respects 



AFi = AEi - AV 



(C12) 



as before. On the other hand, if we further grow a sec- 
ond topological loop so that the torus is encircled both 
directions, at the step where it goes around the torus in 
the second direction we again increase Ei by 2 but Vi by 
1 and Fi by 0. Thus Eq. (C7) remains true also with 
periodic boundary conditions. 
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